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HEAT TRANSFER THROUGH TURBULENT FRICTION LAYERS* 

By H. Reichardt 

SUMMARY 

The "general Prandtl number" Pr ' = Pr, aside from 

A 

the Reynolds number determines the ratio of turbulent to 
molecular heat transfer, and the temperature distribution 
in turbulent friction layers. A^ = exchange coefficient 
for heat; A = exchange coefficient for momentum transfer. 

A formula is derived from the equation defining the 
general Prandtl number which describes the t emperatur e . as 
a function of the velocity, Por fully developed thermal 
boundary layers all questions relating to heat transfer to 
and from incompressible fluids can be treated in a simple 
manner if the ratio of the turbulent shear stress to the 
total stress '^t/T in the layers near the wall is known, 
and if the A^/A can be regarded as independent of the 
distance from the wall. 

The velocity distribution across a flat smooth channel 
and deep into the laminar sublayer was measured for isothermal 
flovr to establish the shear stress ratio '^t/T and to extend 
the universal wall friction, lav/. The values of ^t/r which 
resulted from these measurements can be approximately repre- 
sented by a linear function of the velocity in the lam.inar- 
turbulent transition zone. 

-The effect of the temperature relationship of the mate- 
rial values on the flow near the wall is briefly analyzed. 
It v/as found that the velocity at the laminar boundary (in 
contrast to the thickness of the laminar layer) is approxi- 
mately independent of the tem-oerature distribution. 

The temperature gradient at the wall and the distribution 
of temperature and heat flow in the turbulent friction layers 
were calculated on the basis of the data under equations (s) 
to (4). The derived formulas and the figures reveal the ef- 
fects of the Prandtl number, the Reynolds number, the exchange 
quantities and the temperature relationship of the material 
values, 

*"Die Warmeube rt ragung in turbulenten Reibungsschichten. " 
Z.f.a.M.M., vol. 20. no. 6, Dec. 1940, pp. 297-328. 
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That the form of the v/all and the pressure drop affect 
the results is illustrated by the variation of the thermal 
"behavior of the friction layers in the pipe, channel, and 
flat plate. 

After a discussion of the different definitions of 
the heat transfer coefficient a new formula for the rate 
^of heat transfer is given "based on the maximum temperature 
difference. The new equation differs from that offered "by 
Prandtl by an additional term, that allows for the conditions 
in the laminar-turbulent transition zone. 



IHTHODUCTION 



A survey of the • lit eratur e on heat transfer in tur- 
bulent boundary layers discloses that the problem has been 
treated in numerous studies (reference l) • Because of its 
t echnologica^l importance, the number of experimental projects 
in which empirical or semi-empirical formula,s established for 
various conditions and for various applications preponderate. 

The theoretic-.il principles are but rarely treated. The 
1 i t 3 ratur e t he r ef or e contains only a few general formulae* 
In .G-ermany the formulas by Nusselt and Prandtl are most 
generally utilized. In the English literature it is custom- 
tATiT to introduce the Reynolds analogy which upon general- 
ization by G. I. Taylor loads to approximately the same 
results as the Prandtl theory. 

The -theories to dat.e are based on simplifying assump- 
tions, such as do not usually obtain in reality. The derived 
expressions therefore roqiiired extrapolation based on exper- 
imental results, the extension- extending beyond the original 
range of validity. The practical point of view v/as maintained 
in arranging the s emi-cunpir ical equations and questions of the 
physical significance became secondary. 

The research programs in heat transfer involving many 
technically important special cases in the turbulent region 
fail to allov; the deduction of a general theory v/ithout 
limitations. The solution of this problem is very closely 
related to t he r e s ear ch * o f ' t he flow processes in direct 
proximity of the v/all* 

Before proceeding to an . analy s is of these questions a 
brief survey of the aval labl e t heo r e t i cal contributions should 
beofinterest. 
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REVIEW 0? PHSVIOUS • CONTRIBUTIONS 



The school of Nusselt has made great strides in the 
study of heat transfer prohlems "by the use of tne theory 
of similarity, particularly in arranging the various suh- 
divisions in reasona'ble order. The great technical im- 
portance of the model studies is that it does not require 
the exact knowledge of the individual processes 5 and that 
sim^ole formulae are obtainahle for practical use, even in 
complicated cases. But since no details of the physical 
mechanism are secured the results can "be of a preliminary 
nature only. 

Reynolds (reference S) attempted to define the rules 
of heat transfer from the point-to-point vari \tion of the 
flow pattern. He proceeds from the assumptijn- that the 
turhulent mechanism of heat transfer is the saMe ac the 
mechanism of the momentum transfer. But his considerations 
are still incomplete for practical application and only 
through supplementary considerations by Taylor (references) 
and Stanton (reference 4) v/ere the results of Prandtl 
accomplished, 

Prandtl (reference 5) also starts from the assumption 
that heat and momentum are transferred oy the same mechanism^ 
A complete analogy "bccween these phenomena does not exist, 
however, unless similar "boundary conditions obtain, when 
the nondimons i onal p*Cp/X (termed the- Prandtl number, Pr) 
is equal to unity, and when the pressure drop is ne gl i ..<:^i bl c 
(as, for instance, in flow pa,st a flat plate). In contrast 
the momentum transfer with pressure drop (pipes and channels) 
is described by equations' which differ from those of heat 
transfer and momentum. 

In order to treat the technically important cp-se of flow 
through a pipe Prandtl postulated fictitious heat sources in 
the stream, by means of which a sufficient simila^rity of the 
equations of heat and momentum transfer is obtained. The 
Reynolds concept was taken, that is, that in a very thinlayer 
•near the wall practically all of the transfer is by molecular 
action and that outside of this layer only the turbulent ex- 
cha,nge. mechanism is effective, v/hile the molecular conduc- 
tivity may be neglected. 

The heat source postulate then leads to a simple equatio 
betv/een the heat transfer and the resistance to flov/, which 
may be vrritten in the form: 



4 



HACA Technical Memorandum No. 1047 



where 



o-^j, = — = X — (1) 

^ni 1 ^ ^^(Pr -1) 



a^- heat transfer coefficient referred to the mean 

temperature 

, Cp specific heat 

density of heat flow at the wall 

T shear stress at the v/all 

0 

u mean flow velocity 

m 

Ug^ velocity at the ""boundary" of transition from 

laminar to tur"bulent flow 



To use equation (l) the ratio '^a/'^m ^"^st "be knox^n^ 
In the absence of experimental data of the extremely thin 
wall layer, Prandtl (reference 6) used the following reason- 
ing to evaluate u^^. In the laminar layer a linear velocity 
increase exists, the slope of which is fixed by the shear 
at the boundary. In the turbulent core the l/7th power 
law holds for Reynolds numbers belov/ 10^. The plane in 
which the two velocities coincide is called the boundary 
betv/een the laminar layer and the turbulent zone. The 
exact determination of the boundary velocity u^^ is to 

follow from the heat measurements. 

The heat transfer data available to date indicate 
that Prandtl' s formula do3s not hold for large Prandtl 
numbers^ In consequence there have been proposed various 
corrections to this formula in order to meet the require- 
ments of practice. 

• The basis of the discrepancies lies in the idealiza- 
tion of the transition from laminar to turbulent flow. 
This transition is naturally continuous, hence an inter- 
mediate layer exists in which the viscous and turbulent 
shear stresses are of the sam.e order of magnitude. Since 
the transition to turbulence occurs close to the v/all, it 
has not been possible so far to measure, the velocity dis- 
tribution in the intermediate layer with sufficient accuracy. 
Von Karman (reference 7) has estimated the exchange conditions 
in the transitional region based on an extrapolation of 
iJikuradse^s velocity measurements in the direction of the 
wall. Based on his postulates, von Karman gives the formula 
for the heat transfer coefficient as: 
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— = 1 + [(Pr-l) + t In (l+c (Pr-l))] (s) 

(X - ^ 

where 

ai heat transfer coefficient for Pr-l 

a, "b , c, constants (reference 8) 

An improvement of the theory has been carried. out by 
Taylor (reference 9). Starting from the postulates of 
Reyno Ids-'Tay 1 or , the latter discusses the error of the 
anal ogy • between iieat transfer and momentum transfer for 
flow accompanied by pressure drop* Taylor calculates the 
temperature profile which corresponds to a velocity profile 
measured by Stanton at Pr = 1., The temperature profile 
differs somewhat from the velocity profile, that is, the 
temperature gradient at the wall (and hence o.lso the heat 
transfer coefficient) is lov.rer by several percent than the 
wall velocity gradient. 

Of great practical interest is the variation of the 
heat transfer coefficient for non- i s o t hermal flow in vrhich 
the material values vary v/ith t em.pe ratur e . Apparently this 
problem has not yet been solved analytically. The theories 
to date imply isothermal f 1 ov/ (material properties not a 
function of temperature or space). Since large temperature 
differences do occur in practice, the proper mean magnitudes 
of the material properties are introduced into the isothermal 
expressions ( r ef or enc e 10 ) . 

V/hile the present report was in the press, two further 
articles dealing with turbulent heat transfer have appeo.red, 
one by Mattioli (reference 10a) and the other by Hofmann 
( reference 10b) . 

Starting from special theoretical concepts v/ith respect 
to the turbulent mechanism, Mattioli extrapolates the tur- 
bulent velocity distribution into the semi-lam.inar zone in 
order to deduce from this velocity concejjt the presumption 
of equa,l exc hang.e quant i t i e s for heat transfer and momentum 
transfer the mi^gnitude of the turbulent heat transfer. A 
careful analysis of the difficult derivation shows that the 
important phenomena near the boundary are not adeq^uately 
definr.^d. In addition to the semi-lnminar layer there is 
presumed to exist a v/all layer (v/hich is established from 
the heat transfer mo:-^ surement s of 3\ihne and from other fluid 
flov/ measurements m.entioned above) which is much greater in 



6 



KACA Technical Memorandum No. 1047 



thickness than the laminar layer, Mattioli is therefore 
forced to assume a substantial exchange in his wall layer. 
Since the Mattioli theory connot describe accurately this 
exchange near the wall, the temperature change in the wall 
layer is put proportional to (Pr)^, where m is established 
from heat transfer measurements. 

It is v^orth noting, however, that Mattioli quantitatively 
allo^/s for effect of temperature on the viscosity. For this 
purpose a generalized dista.nce parameter is introduced in a 
manner similar to th.^t employed in the present report (see 
epilation (30) ) . 

Hofmann calculates the temperature distribution and 
the heat transfer coefficient with special consideration 
of the laminar layer whereby the usual simplifying postu- 
lates are retained. The concept of a thermal boundary 
bet'-'een the turbulent cord and a, boundary layer is also 
adopted and the thickness of this layer is discussed. In 
contrast with von Karman, pr ogre s s is made in that the lam- 
inar layer thickness for high Prandtl numbers is introduced. 
The arbitrarily chosen velocity distribution near the wall 
lies above the test points of the present report. 

The position taken by Hofmann that the heat transfer 
depends solely on the velocity distribution and on the 
Prandtl mixing length requires a correction. Basic to 
every theory is a hypothesis of the turbulent diffusion 
of heat. If the ratio of the exchange quantities for heat 
and momentum transfer is chosen (Hofmann tacitly presumes 
the identit^^ of these quantities), then the laws of heat 
transfer follow at once direct from the velocity profile 
without the aid of any turbulence theory, hence without 
the help of the Prandtl mixing length, which in consequence 
drops out again in the course of the Hofmann c alcul a^t i on . 

THE PROBLEMS 



In order to avoid subsequent corrections and to pre- 
sent the hydrodynamic theory of the turbulent heat transfer 
coherently the following assignments are to be solved: 

1. To derive a general equation for heat transfer into 
v/hich the technologically important boundary con- 
ditions and the flov/ phenomena, particularly in 
the transitional layer, co.n be introduced. 
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2. To measuro the flow processos near the wall for 
technologically important casos, particularly- 
smooth surfaces, rough surfaces, alm.ost isothermal 
flov;, non-isothermal flov/, and so forth. 

3« To introduce the obtained data on wall a,djacent 

flow into the general expression to huild special 
. formulae v/hich can he checked hy heat transfer 

moixsur ement s • 

The following statements are made relative to these 
prohlemsi 

She presentation of the general theory should he clear 
from a physical point .of viev/ and it should he simple in 
order that it may be utilized in practice. 

The derivation of a generally applicable equation for 
heat transfer is carried t.hrough in a simple manner. In 
contrast the measurement of the flov/ distribution near the 
wall presents considerable difficulties.. In order to obtain 
practical test data especially thick -boundary layers are 
essential* This requirement implies large flow sections 
and low flov/ velocities, tiiat is, low dynamic pressures 
and low pressure drops must beraeasured* 

The conditions becom.e complicated if the flow is not 
isothermal. Through the influence of the temperature field, 
not only the material properties but also the flov; phenomena 
are changed. 

The presence of roughness introduces further complica** 
tions. It is true that flow on rough v/alls has been exten- 
sively studied and the lav^s of the "nuclear flow" in pipes 
are v/ell known but there is no dependable knowledge of the flow 
processes near the wall between the protuberances. 

The experimental exploration of the flo--/ distribution 
near the wall is a broad field of research which can only be 
accomplished piecemeal. The author first explored the data 
available near the wall. while these studies are not com- 
plete, they have progressed far enough to enable a theoretical 
treatment of the heat transfer at a smooth wall. 

A particularly important sub-task consists in checking 
the applicability of the theoretical formulae by means of 
heat transfer measurements as the theory contains postulates 
relative to the mechanism of heat transfer i^rhich require 
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confirmation "by experience. If necessary the theoretical 
assumptions must "be modified to fit the experimental facts. 
The heat transfer measurements can be employed with great 
benefit to clarifj'- the questions of turbulence structure* 



THE PRANDTL NUMBER 



The hydrodynamic equation for the continuity of heat 
flov/ (equation (47)) is not sufficient for predicting the 
temperature distribution in the friction layers. It re- 
quires another equation for the temperature which takes 
into account the requirements of the system under consider- 
ation* (This temperature equat ion, 1 ooked for, places the 
continuity of heat flow equation (equation (47)) in the 
position of a special condition that must always be satis- 
fied.). 



The Prandtl number Pr 



governs the form of the 



temperature profile* It is logical therefor e , to begin 
•v/ith'the Prandtl number concept. To secure a *di f f e rent ial 
equation necessitates a determinating equat ion for Pr that 
holds for each point in the fluid. Since the individual 
factors in Pr have "point" significance, the derivation 
of such an equation is possible. 



Let q equal the density of the heat flov;, and T 
the shearing stress of the density of momentum transfer. 
Assume that the heat flow and momentum flow act in the same 
direction {±y) at a given point, which is perpendicular 
to the mean velocity u (time average) at this 'point. . In 
the system under consideration y is measured perpendic- 
"ular to the wall- and u parallel to it. 



The total momentum t 
the molecular transfer and 

exchange motion. The same 

Hence 



consists of a portion by 
a portion by the turbulent 

holds true for the heat flow* 



T = 



m ^ t 



(3) 
(4) 



vath 



coefficient of viscosity 
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}\ thermal conductivity . . . . 

A' exchange coefficient for momentum 
A^ corresponding coefficient for heat 
Cp specific heat 

u, T time averages of velocity and temperature, respectively 

u', v' velocity f luc tuat ions in the x, y, directions, respec- 
t i ve ly 

T ' corresponding t erap.;; r atur e flue tuat ion 

'^m = ~ (5) 
dy 



du 



. = A — ' = - o u ^ V ' • 

^ dy .' 



(6) 



^^m 



= ^ — (?) 
dy ^ 



^^t = ^p — = - P 

- dy • 

The coefficients p.,. A, A, A^, are defined hy these equa- 
tions. There presentation of ^ and - in terms of the 

fluctuating components is for the present irrelevant, "but 
Will be clarified in Chapter 9. 

Equations (5). to (8) then yield the . f ol lov/ing proporti 
^■t ^t - 

— = Pr — = Pr » — . . . (9) 

Accordingly the ratio of turbulent to molecular heat flow is 
proportional to the ratio of- turbulent to molecular shear 
stress. The proportionality factor is • Pr » = (AQ/A)Pr, a 
quotient which is called "general Prandtl number." 



on 
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Equation (9) thus leads to an extension of the concept 
of the Prandtl number for turbulent flovr with P* = A^/A Pr 

instead of Pr. Only in the case v/here the exchange coeffi- 
cients are the sane for momentum transfer and for heat 
transfer will Pr and Pr ' be equal. ' 

Since equation (9) refers to flow in which turbulent 
and molecular shearing stresses act, it is particularly 
'suitable for the representation of the physical phenomena 
in the transitional layer. The treatment of the heat 
transfer in the present report therefore starts from the 
transitional flov/, the "fully turbulent" core and the 
laminar motion at the wall being treated as special cases. 
(In the proximity :of tha v/all the exchange mechanism per- 
pendicular to the v/all is not possible; therefore the 
turbulent friction disappears and the momentum transfer is 
accomplished by internal friction only. Because of the 
turbulent pressure fluctuations, the stream velocity near 
the wall also experiences f luctnat ions • The continuity 
of this fluid flow is largely maintained by the lateral 
transverse fluctuations, so that the wall flov7 glides 
practically parallel to the surface. In this sense the 
viscous wall flov/ is "laminar.") 

A picture of the physical significance of the Prandtl 
number is best obtained by observation of the transitional 
layer for extremely high values of Pr ' (very viscous fluids) 
In this case practically only turbulent heat transfer exists 
(q^» q^) at those places in the transition region where only 

sm.all turbulence exists ( '^x ^< ). In this extreme case the 

^ m . 

molecular heat transfer is so small, that even a slight con- 
vection signifies a form of "short circuit" for the heat flow 
Therefore the temperature profiles for high Prandtl iJumbers 
are " smoothed . " 

Even for the special case of Pr ^ 1 and q/T = constant 
the temperature profile can be fixed readily. It is 

'^t- q ^m ^ dT ^ 
— = — - or — = "IT" = • In this case theprofiles of the 

^m m m 

t e:a-peratur e and velocity agree with each other. (The con- 
dition that q = T is v;ell satisfied in the friction layer 
of a flat plat e . ) 

THE GENERAL TEMPEHATUHE EQUATION 

The temperature distribution follov/s from the equation 

of the molecular heat stream qm = A qf must, therefore 

dy ^ 
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be replaced by ( q. - q^^^) in equation (9). Introducing the 
ratio q/T equations (5) and (?) then gives: 



i-m 



X ^T q/r 



du 1 + (Pr« - l)'^t/'^ 



(10) 



The boundary condition at the wall is to be introduced in 
this general equation. That is, 



(11) 



0 o 



and v/hen augmented by (H), equation (lO) is integrated t< 



The factor ( — j is determined by extension of the integral 
\du/o 

over the total velocity field of the friction layer. 



/ dT^^ /^dT 

The temperature-velocity quotient — = ( — - 
^ ^ \du/o Vdy/ 



^duN, 
dy^o 

is a measure for the amount of heat transferred to the v/all. 
The heat transfer at the v/all is obtained from the temperatur 
dis tr ibut ion. 



Equation ( 1 2 ) , al t hough designed to calculate the temper 
aturo distribution, has ge^neral application. The considera- 
tions so far are based solely on knov/n definitions properly 
rearranged and combined and no special assumptions relative 
to the flov; have entered the computations except the bound- 
ary condition of a laminar wall layer. 

The above derivation indicates that a general result can 
bo secured without employing the hydrodynamic equations 
(46) and (4?)). This is due to the f3vCt that the basis o 
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each theory, entirely independent of the method of calcu- 
lation, is a postulate related to the exchange mechanism. 
(For instance T' is set proportional to u', or more 
generally/" A^ = A, or as in the case in point, Aq/A is 

to "be determined later^) This simple hypothetical content 
of the theory is seen also from the presentation of the 
simple Prandtl analogy for the present subject is treated 
in such a manner as to make this step possible. 

To coijiplete the tem.perature equation (12) the magni- 
tudes A^/Af q.T^/q^T and t"^/ t must be known. (These 

quotients are introduced later in order that the effect of 
each postulate may be observed independently. Also the 
various deviations between the theory and experiment reveal 
at a glance the direction in which the assumptions must be 
modified. ) 

In order to carry through a calculation Aq/A is 
assummod to be constant. The value of Aq/A is to be 
determined from experimental data. 

The quotient qT^/q^T cannot be fixed arbitrarily. 

The heat stream- q • • is -relat ed to the temperature T 
through the differential continuity equation (see equation 
(47 ) gonerali z ed Pour ier-Poi s son) of hea.t flov/. But a first 
approximation of the temperature distribution can be obtained 
by assuming that the layer for heat tr^insfer is of about the 
same thickness as the friction layer. 

In this case the heat flow disappears v;here the shear 
stress is zero, while on the v;all ^/ = 1 '^/'^o ~ "^^ 

Thus- the total range of the friction layer can be expressed 
V/ i t h 



^ = 1 + k (13) 

... 

where k is small compared to 1 at least in proximity of 
the wall. • 

In the turbulent friction layers the velocity gradient 
is steep near the wall. The largest part of the velocity 
region u lies in a zone wherv3 k is small. So for the 
inte^-ration T over u of equation (is), (^"'^o/^o'^^ ^ 
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may "be .put in first approximation. (In the entrance zone's 
where the wall, t'emperature changes suddenly this approxi-" 
mation is not pb s s-ihle/; ."H'or such cases the heat boundary 
layer is much thinner than the friction layer and. it there- 
fore plays an • important role in the variation in heat flov/. 
Thermal entrance le.r.gths in e.xisting friction layers are 
quite short hov/ever (see Latzko, Zoa.M.M., Bd . 1 (l92l) 
p. 268) , so that when assuming (j-i A q/^o^ ) 1 and, 
T^/t is known, the integration can he completed. 

This procedure yields a first approximation of the 
temperature distribution hy means of which the heat flow 
can ho evaluated. The heat flow distrihution then affords a 
second approximation for the temperature di s tr ihut ion which 
is practically adequate for the case of constant material 
propert ies . 

Several quantitative conclusions can "be drav/n from 
equation (l2) relative to the tempv^rature profile of various 
friction layers which coincide approximately with the stress 
quotient ( / '^) (such as, for example, the flow through a 
pipe, channel or flat plate) at equal Reyno Ids. number s , 
v/hore the velocity distribution obeys tho "universal law"). 

At the flat plate ( ) = 0 and likewise ( — ) =0. 

Vdy/o \cLy/o 

The assumption k 0 is therefore vrell satisfied over the 
greater p.'-rt of the velocity field of the flat plate. l^o 
appreciable differences obtained here between the first and 
second approximation and the final solution of T. (Sven 
though T and q are very similar at the plate, they are 
not coincident, for q .depends' on Pr' vrhile t does not. 
Therefore there will exist for the plate, a small difference 
between the actual temperature profile and the first approxi- 
mation of T.) 

Por flo'.7 wi th -pr es sure drop a far f rom . negl igible 
difference exists betv/een the second and first approximations 
(that is, between the actual profile and that of the "plate 
profile" of the t erap^ ratur e . ) By pressure drop {d'^ / d^y) ^ < 0 , 

but at the flr^t channel v/all (dq/dy)^ 0 and in the pipe 
(dq/dy)Q> 0. Hence it follows* for the temporature distri- 
butions that the pipe profile differs more from the plate 
profile than the channel pr o f 1 1 a , • t hat is, according; to equa- 
tion (l:^), the temperature rise at the pipe wall is flatter 
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than at the' channel v/all and even inore so than at' the fla^t. 
plate '(sec fig. 4). Por a quantitative treatment of model 
problems it is advisable to integrate the temperature 
'-equations by sections, that is, the laminar section, the 
.trans it iano.l region and the real turbulent layer. The bound-* 
a,ry at the end o'l the. laminar zone is designated by the sub- 
script a, and t.ixe .beginning of the turbulent layer by b^, 

whence, after introducing the substitution equation (is),' 
equation (l2) gives . . . ■ .."^ 

^ /du^ ^ r ^-^^ ^ ■(Q< T<TJ (12a) 



VdT'^ 



rs ^ U A 

0 0 ^0 



/. N % n ^(1+ k) 

/duN r r . / s 

0 a T 



(^V = r.r. ^ ^ ^^^^ . ''^ • (T,,<T<9) (120) 



t. 



'(9= ma.x, temp* difference between the wall and the f levy- 
ing fluid.) In addition it should be observed that k may 
be disregarded for the lam inar • r eg i on • A general disregard 
of k in the main fluid stream is not tenable.- The sub- 
script t indicates a mean value for the turbulent region 
( f oriaed over u) , . 

In the actual turbulent region it is to be noted that 
for small values of Reynolds numbers '^t/^ considerably 
sm^.ller than unit (see fig. 2), correspondingly '^m/'^ 

not negligible. However to an apior oximat ion ('^t/'^^ ~ 
('^t/"^t ~ constant. The point where '^^/'^ = ^'^t/'^^t 
the turbulent boundary designated with b^. 
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To utilize the temperature equations the variation 
of (t-^/t) and the "boundary velocities u^ , ut) .(and corre- 
spondingly u-u ) must he knov;n; i^hi s involves the flow 

0 

distribution noar the v'alls, v/ith limitation to the processes 
at the smooth v/all and tjD j'lov/s obeying the ..uiU.Y.ersa.l_veloc i ty 
d i s t r ihut ion e quat i on • 

VELOCITY DISTRIBUTION AT A SMOOTH V/ALL 

The measurements by -iJikurad-se (reference ll) have shown 
that the turbulent velocity distribution can be approximately 
represented "by t he ' following equation; 



— = 5.75 In y* + B (14) 
u* 

v/here the dimens ionlo ss she ar in r^' s t r e s s velocity is defined 
Dy 




(15) 



and the d imen s i o n 1 6 s s v/all distance by 



y* = -ir (16) 

The constant B depends on conditions at the '"/all, 
Por smooth walls B is approximately 5,5. -Equation (14) 
is a straight line "on s em i- lo gar i t hra i c pape r as shown in 
figure 1. 

The velocity distrioution for the laminar wall layer 
can £ilso be represented by means of u/u*** and y*. He- 
arrangement of the Poiseuille equation results in 

-ii- = y* f 1 - (17) 
u* V 2 



where n = y/ r and r 
for the equation in the 



is' the radius of the pipe or channel 
sub-layer. In general the laminar 
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layer is so thin that ti/2 may be neglected compared to 
unity and therefore practically 



(17a) 



Equation (l7a) is therefore the universal equation 
for the velocity distribution in the laminar zone. It is 
shovrn in figure 1 as the curve which passes through the 
point logy* = 1, u/u* = 10. 

The flow conditions in the transitional layer are not 
very well established experimentally. This sublayer 
adjacent to the wall is usually so thin that accurate measure- 
ments of the velocity can hardly be made. The closest v/all 
proxim.ity was probably reached by Stanton with his surface 
tube (reference 12), But even these test data are insuffi- 
cient for the present arguments. 

Since the application of our theory is predicted on 
the knowledge of the shear stress ratio (t^/t) in wall 
proximity^ a wall layer of such thickness was required 
as to render a measurement of the wall flow possible. 

The thickness of the laminar layer yg^ and the bound- 
ary velocity Ua are fixed by definite values of u/u*** 

and y*; y^ increases v/ith decreasing u* according to 

equat ion (16) ♦ The reduction of is limited by the fact 

that at too low shear forces the critical Reynolds number 
is undercut and so the entire flow becomes laminar. It is 
therefore appropriate to introduce the Reynolds number 
'. /Re = "^m^/^ in the place o.f u*. Then the thickness of 
the v; a 11 layer is 



y 



a 




(18) 



where the so-called re,sistance coefficient i isdefined in 



the usup.l manner as J 




(^ decreases slightly with Re). 



(19) 



NACA .Technical Memorandum Ho. .1047 



17 



T-he t-hickness of the wall layer .^rovrs v/ith the 
diameter d of the pipe or cha.nnel and decreases v/ith 
,the P.eynclds- number. For a given Reynolds numhor 7q,I^ 
is independent of the choice of flowing mec^ium. 

The important numher 7 y the exact value of which 
...is not yet knov/n, lies "below 10 according to available 
measurements. The critical Reynolds number is 3000, and 
the corresponding t ^ 0.04. Herev/ith 

This equation reveals that even for the lovrest possible 
Reynolds number the stream diameter must, be fairly great 
in order that a probe can be introduced into the laminar 
wall layer. (With considerations to the influence ,pf the 
wall in the probe, the wall layer should be as thick as 
possible.) 

But the achievement of a sufficient boundary layer 
thickness by increasing the stream diameter introduces 
fundamental difficulties. If the increase in diameter is 
to achieve the purpose desired, the Reynolds number miay 
not be increased (see equation 18). This means that the 
velocity must be decreased in the same })roportion as the 
diameter is increased. As a result the dynamic pressure 
and the pressure drop are reduced quadrat ically with the 
stream diameter, that is, v/ith the boundary layer thickness 

This is exemplified at the dynamic pressure of the 
mean velocity Uju, for which introducing Re = du^/ V, we 
get 



^o'fix the order of magnitude of this dynamic pressure 
several numerical values are inserted. Let Re ^ 3000, 
d = 25 centimeter .(utilized in measurements reported 
later in this pajjer). For air as the flowing medium 
(m. = 1.8 X lO-*", p = 1.2 X 10-^) it is 

^ - 2 X 10-'^ mm. H^O 



2' o d^ 



(21) 
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The yjressure drop to be measured is of the same order of 
magnitudes and is equal to p Ujjj^ at Re = 3000 for 
Al = 100 d, Por flovjing water th e s e • pr e s sur e s ar e about 
four times larger and only when utilizing a viscous oil 
does the magnitude become equal to 1 millimeter of water* 
If these pressures are to.be measured to within 1 percent 
then the sensitivity of the manometer must be in the range 
lu~^ to lO'"^ millimeters of water. 

The problem of precise measurement of the flo\^ phenomena 
close to the wall for non-isothermal flow involves the tech- 
nical difficulty of measuring extremely small pressure differ- 
ences. (The velocities can be determined without the use of 
pressure measurements. In the boundary layer itself a hot 
wire anemometer or a thread anemometer can be used in place 
of a pitot. tube. These devices. must be calibrated and the 
calibration at best depends on pre s sure, measur ing devices* 
In addition pressure drop me asur e-ment s are desired to check 
the effective shear s tr e s s • ) Such measurements can be made 
with the micromanome tar designed by the author v/hich has an 
upper limit of sensitivity of 10"" millimeters of water 
( reference 13) • . 

The turbulent flov; measur emsnt s reported here were made 
in a rectangular channel 35 cent imeters . high , 1 m.illimeter 
v/ide .'^.nd 16 mill imeters .long and with a maximum velocity of 
80 centimeters per second. jine pitot tmbes and hot wires 
v/ero utilized. The hot v/iro anemometer was ca^librated in 
the parabolic distribution of a 3 centimeter high x 30 centi- 
meter wide laminar channel in v.-hich at similar distances from 
the Weill, the same Tq obtained c'ls in the turbule.nt channel* 

was evaluated from the maximum velocity as well. as from 
the pressure drop. 

The measurements vrere made very difficult because the 
lov.r velocities were easily disturbed by external causes. 
For instance, small teranerature diffprences between the 
air stream and the v/all (induced by unavoidable fluctuations 
in room t em;pe ratur e ) caused observable changes in the veloc- 
ity distribution. Therefore the turbulent velocity profile 
was almost alvra;'''s slightly unsymme t r ical and thus u v/a^s 
different on the upper and lov;er wall. 

On the top of t h^.t . t he r ecord'ed pressure drops yielded 
an average u* v;hich was too great as compared with the 
results of other authors. The channel flow obviously was 
not completely developed (in a tube the length of 64 diam- 
eters v;ould have been sufficient.) But since the pressure 
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drop was held constant for all test points, it was 
possible to determine the raean u* "by comparison with 
indisputable measurements of other authors at higher 
values of y*. Por this purpose measurements of Nik*- 
uradse v/ere utilized omitting those for v/hich the wall 
correction was questionable,- (Similarly the measurements 
of the Stockholm report which do not lie in the range of 
others of Nikuradse^s measurements and are obviously too 
high, have been om.itted.) 

The results of these m.easur ement s necar the boundary 
are shov/n in figure 1. The u/u* points approach the 
laminar curve very gradually. It is reached at approxi- 
mately u/u* = 1.5, a value v/hich is substantially lower 

than that usually assessed. 

The value Ug^/u* =1.5 is however still uncertain 
and it must finally be based on much more accurate measure-^ 
m. e n t s • 

It is also true, that an accurate determination of 
the limit where du/dy (du/dy)^ is not possible from 

velocity measurements. For this jpurpose heat transfer 
measurements at high Prandtl num.bers will serve better 
to ci.etermine the l^-^ininar boundary. From heat transfer 
measurements by Buhnc- it would appear that u^/u* is 

somev.^hat larger than 2. 

The recorded velocity distribution in the transitional 
region can be approximated at: 




(22) 



b-a 



b-a 



wher e 



a = 



u* 




and 



n,^ is the velocity at the laminar boundary and u-^ a 
suitably selected velocit-- at the turbulence boundary. 
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As "bofore r\/2 can be neglected compared to 1 and 
equation {22) then becomes a universal law. In figure 1 
equation (22) ic presented for a = 2, b = 15. The curve 
is. dashed above u/u* 15, where it loses its physical 
significance. The measurements are satisfactorily rep- 
resented by this equation. • 

It remains to be explained why the velocity distri- 
bution in the transitional zone was appr oximat ed " by 
equation (22) although som.e other similar function had 
been possible. 

The ratio t^/t is required. To fix the ratio, 
differentiate equation (22): 

/ \ 

d ^ ^ 1 

VW_ 

(l - ri) dy* b - a 



wnere 

d ■ T 

- ^ ^ -Ja (24) 



dy* qvl"^^ dy 



0 



The total shear stress for developed flow v/ith pressure 
d r 0 p i s ■" ' 



T = ro (1 - T]) (25) 

(near the wall one may set ~ '^c''* Then, solving, 
one obtains : 



' m b - u 



^b ^a 



u - u. 



^b - ^a 



(26) 



(2?) 



Since eouation (22) is confirmed quite well by the measure- 
ments for '^^< ^< the true variation of should 
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not differ much from equation (27). An uncertaintly 
exists, of course, at the limits^ a and .h, 

There- fol-lows from e-iji^uation (27) "by introducin;^ the 
often used ratio = u/U 



JJ 

t u' 



(27a) 



where U is the maximtim velocity at the edge of the friction 
layer* U/u* is, in accordance with ec^uation (l4) a function 
of r* (that is, the value of y* at the "border of the 
friction layer "based on the distance r- from the v^alk). 
The relation between Re and r* is given "by the identity* 



U U . U r u* 
R = r* = 

wher e 

Ujp^ mean velocity 

— = a function of He 

U a 

Figure 2 shows '^t/'^ different values of QO 

v/ith Reynolds numbers £ts the parameter as calculated from 
e(]^uation (27a) for the transitional layer and by equation 
(14) for the turbulent region. The constants a and b 
were chosen at 2 and 15. 5,. respect ively. 

The actiial '^t/'^ distribution no doubt differ:S from 
that shown in figure 2 for small values of ^. But the 
difference "between the velocity di s tr i'but ion as expressed 
"by equation (22) and the lamina.r curve is less than the 
scatter of the experimental points (see fig. l) so that 
nothing certain may be said relative to the actual ''"t/'^ 
variation near the laminar "boundary nor of the laminar 
boundary itself. 



(28) 
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Details of the -variation of '^t/'^ P^ay at first no 
part. In contrast v/ith earlier vrork in which the friction 
layer was dividod into tw:o regions in which '^t/'^ varied 

from 0 to 1, it should for the first suffice to approxi- 
mately describe the processes in the transitional zone. 

The earlier division of the friction layer into a 
laminar and a turhulent region is indicated "by tv/o vertical^ 
lines in figure 2. The dotted line represents Re = 4 x 10 
and a = 2 and the dot dashed line represents Re = 4 x 10^ 
"but a = 8.8, which is the value chosen by Prandtl in 1928. 
t^/t was defined as zero up to a = 2 (or 8,8) and unity 
for greater values of y*.) At high Prandtl numbers where 
the transitional layer can be regarded as pa,rt of the tur- 
bulent zone core with respect to heat transfer (see equation 
(9)) a = 2 is in good agreement while a = 8.8 results in 
a heat transfer rate v/hich is too small, 

THE EFFECT OF TEMPERATURE RELATIONSHIP OP THE 
MATERIAL VALUES ON THE FLOW PHENOMENA 

If the material values are functions of the temperature 
then the flow distribution across a section will be changed 
as mentioned above and also reduction in temperatures in the 
direction of flow causes hydrodynamic changes for all fluids 
v/hich are compressible. In this instance,, in principle at 
least, there exist no velocity profiles which are similar^ 
the same statement holds for the temperature profiles both 
cons ider ed as a function of length. 

Since the magnitudes of the temperature differences and 
the differences in the temperature coefficients of each prop- 
erty enter into the evaluation of the profiles, a general 
solution of the problem is hardly possible and the study re-? 
stricted to the simple case of similar temperature profiles 
which are practically achieved at relatively low differences. 

If the viscosity of an isothermal friction layer is 
changed from to "^s and if the remainder of the vari- 

ables, particularly u* do not change, then equation (14) 
reveals a parallel displacement of the turbulent velocity 
profile (see fig, 3a) with a velocity difference of 

Au = 5.75 u* In -— (29) 

1 



KACA Technical- Memorandum ITo. 104.7 



23 



From this it follows that the viscosity has practically 
no influence in the fully turbulent region, "but affects 
solely the "boundary velocity near the wall,. (The tur- 
bulent flow' slides" at' the wall at a higher or lower 
velocity ec[ual in magnitude to A u,) 

If the viscosity in the turbulent core of the iso- 
thermal flow plays no part its influence for non- i s o thermal 
flow is limited to the effect due to its variation. Vis- 
cosity variations in the turbulent core are not ^^reat, for 
the temperature variations are not great. We may,, there- 
fore, generalize the laws established for isothermal flow 
by omitting the effect of viscosity in the turbulent core 
and by replacing the isothermal viscosity ^ in equation 
(14) by a suitably defined laminar layer viscosity V , 

Recently the resistance measurements, of Hohonczi 
(reference 14) for non- i so thermal flow of hot water 
being cooled in a tube were published. The measurements 
could not be adequately correlated if the f r i c t i on f ac t or s 
were plotted against Reynolds numbers in which the viscosity 
is evaluated at the mean fluid temperature. In contrast 
the correlation is satisfactory if the viscosity v ^ is 
evaluated at the v/all temperature. 

The rest of the discrepancies can be eliminated if 
a viscosity slightly less than that corresponding to the 
wall is employed in the Reynolds numbers. As far as the 
author could determine the results of Hohonczi can be 
satisfactorily correlated and are in agreement with those 
of Bias ius-Nikuradse if the Reynolds number is referred 
to the mean laminar layer \^iscosity 1^1^ and a * is put 

equal to 2. (Hohonczi chose v ^, as the correct viscosity 

due to an error in conclusion f r om s im i lar i ty reasoning 
in which the differential equation for isothermal flow 
was applied to n on- i s 0 1 he rmal flow. In addition the 

values of Hohonczi do not achieve coincidence of the iso- 
thermal and non-isothermal results. Up to this time the 
thickness of the laminar sub-layer was chosen too thick, 
resulting in a sublayer temperature which was too high 
and a value of vl which is too low. Only for one set of 
data at high Reynolds number "v;ill the results yield to 
adequate ccrr elat ion. ) 

Consider next the influence of a uniform viscosity 
variation on the flcv; conditions near tiie wall. If the 
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friction velocity u* is not changed* then the boundary 
velocities u^^ = a. u* u-i^ = h u* are maintained since 

a - and "b are universal constants." That is, only the 
layer thickne s.s es change to 

ya* ^ 2''b* 

= — — ' yi) = — — 



In figure 3a the velocity profiles for unfform changes 
of viscosity are shewn in which, for the sake of simplicity, 
the transitional layer is included with the turbulent core. 
Curve 1 is the original profile. Heducing the viscosity 
yields profile 2 v/ith one^-half the laminar layer thickness* 
Increasing the viscosity by 50 percent yields profile 3 
with a cor responding laminar-layer thickness of 1.5 of 
the original layeri 

The ca.se of a locally vp.r iablo visco s i ty such as obtains 
in non-isothermal flov/ stipulates a generalization of the 
dimensionless distance: y* = yu*/'"D For th3 viscous wall 
layer the follov/ing simple possibility presents itself: 




(30) 



' Tha p.pplicabi lit j'^ of this concept must be established by 
•experiment. But it may be stated that this concept (equa- 
tion (so)) is more satisfactory th^in" the original and that 
one Can predict well those cases in which the property- 
temperature quotient is not too great by. employing equa- 
tion (30) . 

The ratio u/u*' can be genera^lized by re-arrangement 
of equation (3) for the laminar wall layer. 

' p ' d y p d u . 



The following indentity holds for the. laminar layer: 
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-i- [ p d u = f iiZ = (31) 

' ^ 0 ^ 



So yu*/v is replaced "by u* = y*, then u/u 



1 r ^0 



must "be replaced by / o du in' order to preserve 

' 0 0 

the universal representation of the Poiseuille law. 

Under the postulate that there exists a certain •• 
critical number in this r epresonta.t ion , the laminar 
boundar^^ of t he non-i go thermal flow must be at the same 
value as for the isothermal flow ;/^* ~ a. • iTor the ratio 
u^/u* equation (5l) then gives 

ua^^jap^ - ^^^^ 

^* Po+Pa 

if in the first approximation p(T) and T(u) are linear. 
At constant density u^/u* = a as v/as the case for iso- 
thermal flow.- 

in incompressible fluids the laminar boundary velocity 
therefore alway:.s .approximates to the same value a,j^* no 

matter v/hat the viscosity variation in the laminar layer 
may be. The integration limit u^ in equations (l2a) to 
(l2c) can therefore be retained for non- i so thermal flow 
also* 

The dffect of the viscosity?- expresses Itself in the 
thickness of the laminar layer, accordin,?* to equations 
"(30) ^nd (31) : 

a 1 
• y. = ~ f ^' dy* = ^ r ^ d. (^) (33) 

0 0 

These ratios are expressed qualitatively in figure 
3b, that is, for specified values of u* and jj^o st the 
wall. The viijcosity in the sub-layer is smaller in Profile 
2, greater in Profile 3, than the viscosity m-q ^^^^ iso- 
thermal flov; illustrated in velocity Profile 1. 
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It is seen that the parallel shift of the turbulent 
profile in the non- i so thermal case is much less than by 
the uniform variation of The velocity change may be 

expressed approximately from equation (l4) as: 




0 



(34) 



Because of the neglect of the transitional layer this 
shift is less than the true Au, An improvement is 



To further evaluate equations (l2a) and (l2c) 
[iAq/p,qX is substituted for CpQ Pr/cp Pr^ and the 

dimens ionless ratios ^= t^/U, = T/ 6 introduced, 

Tq simplify the calculation a constant Prandtl 
number Pr (-t^j) and a constant specific heat Cp^ is 

introduced for the transitional layer. For the turbulent 
re5:ion itself Cp is equal to Cp^ , a constant* Further 

the ratio of the exchange quant ties is assumed identical* 
(For completely laminar flow the ratio • A^/ A loses of 

course its significance.) 

Defining a mean Prandtl number for the laminar layer as 



1 



possibly obtainable with the integral 



0 




TEE T5MPERATUHE DISTRIBUTION AND THE TEMPERATURE 



GRADIENT AT A SMOOTH WALL 




(35) 



and putting 




T 



(36) 



xTACA Technical Memorandum No. 1047 



27 



it fellows from equations (l2a) to (l2c) upon the intro- 
duction of equation (27a) and the application of ' q ~ 
A^/A Pro and Pr»^ = A^/ A Pr^ that: 

/ deo\, 

^'^'o ^ = ^^^^'a ^ - ^^'l ^a (0 < ^ < (37a). 

V , V c CD / CD-CO 

°V<i^4 V cp- 1-1/Pr'^ V O^-c^^ 

Fi-' — :• = Pr'iCD + e + — r- In - ^ - + :r- (^d-cp^) 

" Vd*/5 ^ a c- l-l/?r'- ^ Cpt ^ 

Jrom equation (37c) the temperature gradient at the v/all 



1 s 



^t ^ 



0 

'DO 



where ei is the value of e for cp = 1 multiplied "by 
^Pt/^po' (See equ-^ttion (36).) 

The quantities CC'„ = a u*/U and ^-^ = bu'^U are 
known functions of r* and He (see equations (l4) and 
(28)) if ^ and b are fixed. On the other hand, the 
material values Pr-j, ■•nd ^^-q; must oe defined more 

ezcplicitlv. (To "be discussed latter.) 

In equations (57c) and (3S) ^--^^ was replaced "by ''^-^ 

because a slight variation v/hile defining the turbulence 
boundary has practically no effect on the calculation of 
the total temperature. (At the boundary position cd^q 

the temp, difference over the turbulent region is greater by 
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^b^^'^lDO f d^A 

, \dQoJo ^^'^^ ^"b» v/hile the .temp, difference 

in the transitional zone is reduced by approx* the same 
amount . ) 

In addition to the very small error terms, the 
expression 



rp - ^^^,) (i-i/Prtj 



t 



is also emmitted, since it is smaller than 0.02 (see fig, s) 
oven for the largest ('^m/'^^t small Reynolds numbers^ 

The number ei accounts for the effect of the variation 
of q/r on the temperature gradient at the wall. Since 
this term is less than unity (see fig. lO) it plays no 
important part except at low Prandtl numbers as is seen 
from equation (38), The integration of equation (36) 
between b and unity usually suffices to calculate ei 
and by this operation Fr ' and Cp disappear: 



e 1 



1 



d CP (39a) 



while the errors due to the inaccuracy of the speci- 
fication of '^t/'^ "^-^ means of equation (27) tend to dis- 
appear for high Pr, for Pr ' = 1 the only term which 
is in error is ei (for '^t/'^ eliminated) to the 

extent that the errors due^ to the material values- can 
be discounted. 

For Pr^ = 1 




(o9b) 



Here the integration from 0 to b for pipe flow is 
appropriate, if the Reynolds niimber is small (in this 
event k in the traiis itionial zone cannot be neglected.) 
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!?or constant ma.terial values D-nd for cd -j^ = ^-^ ^9 

q/T a /t (e 0 ) oqr.ations (38a) to (37c) and (38), 

' "0^ 0 

give the Prandtl approximations: 



a 



/cicp\ 1 , , 

* Wo = + ^ ^^a) (^^ < < 1) (40b) 

Fr' ^ =14- '^O^ (Pr' - 1) > (41) 

(Prpundtl employed Pr instead of Pr ^ , Ho also used 
Uttj and as reference values rp.ther thcan U and 6.) 

The fourth term of equation (36) v/hich accounts for 
the conditions in the transitional layer is particularly 
important for average values of the Prandtl numher. But 
at high values of Pr * the. fourth term is small com« 
pared v/ith the third and for constant material values 
Prandtl equation (4l) is appr ox iniat o ly o"btained again. 

At Pr ' 1, and constant material values equation 
(38) simplifies to: 



IPr^ = I ^ e[ ^^^^^ iVr^-l) (42) 

In flov/ with pressure drop, consideration of the heat 
flov/ distribution which enters into the ei term yields 

a smaller temperature increase at the wall than "by the 
assumption q/ t = q / t . Por instance (d-^/dcp) is not 

equal to unity at Pr * = 1 and He = 4.10 hut in a 
channel is only ahout 0.94 and in the pipe a,ppr oximat ely 
0.91 (see fig. lO). 

In figure 4 (d'^/d(p)^ is pr e sent ed - as a function of 

Pr ' for constant material values with He as the param- 
eter, a = 2 and h = 15. 5^ The solid curves refer to 
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the flow through a pipe* For Re = 4 X 10^ the dot-dashed 
curve is that of a channel v;hile the dashed curves refer to 
a flat plate. 

The asymptotic limit value of (d^'^/d''^) for extremely 
high Prandtl numbers is l/^a* temperature gradient at 

the v/all (dn^/dy)^ is therefore at the most l/^a ^in^^s 
greater than the corresponding velocity gradient, 

loT Pr = 0.72 (air and other gases at room temp.), 
(d^/dco)^ ^ 0.8 at the plate, if A = A (see equation (42) 
and fig. 4). 

Elias (reference 15) has o s t a^hlished , for the flow 
along a heated plate, that the temperature and velocity 
profiles are simil^.r, that is, that (d'^/d'^)^ 1. 

This value for (d^/dcp)^ holds, however, for 
(Ag^/A) Pr 1. Fro'm this it follows that A^/A ^ 1/0. 72 - 1* 4. 

A similar result vras obtained by Lorenz and Priedrichs 
(reference 16) in their experiments v/ith air flowing through 
heated pipe. Re ^ 10^, (d-^/dcp)^ - 0.97. This value lies 

at Pr* - 1.08 as may be seen from figure 4 (equation (l2)). 
Prom this it follows that A^/A = 1.5. 

The question regarding the ratio of the exchange 
quantities, however, cannot be cons idered "solved , hence 
no specified value -^q/-^ v;ill be ascertained* 

Figure 5 illustrates the temperature distribution 
-^(V) for various Pr ' at Re = 4.10"^. The solid curves 
indicate th.e" second approximation for pipe flow, the 
dashed curves represent the first approximations which 
approximately correspond to the temperature distribution 
along the plate. The division into three flow regions 
is indicated by the lines ~ constant, = constant 

(that is, for a = 2.0 and b = 15.5). 

In figure 6 the temperature distributions of figure 5 
are plotted against the d im en s i onle s s v/all distance ri. 
For purposes of clarity only the case of Pr ' = 1 for the 
first approximation (also a near appr^ox. for the plate) 
is presented. This curve also represents the velocity 
profile for Re = 4.10"^ for the universal velocity dis- 
tribution curve approximately holds for the pipe as v/ell 
asfortheplate. 

As to the non-isothermal problem the least trouble 
is in the choice of the value 01 Cp since the specific 
heat varies but slightly with the temperature. In many 
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cases 021:6 may v/rite CpQ = ^pu ~ ^pt ivhich event 

eq.uation (38) is greatly s.implifiea. 

Per a more accurate analysis the approximate range 
of the pertinent temperatures must be knov/n. Employing 
the suhscriptsut ilized to describe the material values 
and arranging the temperatures in the order of increasing 
temperature results in 

I ""a u D u t 

Here, in addition to the material value temperatures, 
the boundary temperatures and as well as the 

mean temperature of the flov/ing fluid x^^J^ are introduced^ 
(The definitions of -T^ and T^ depend on the variations 
of the material values and are very complicated, (see 
derivation of equation (38) But it is not necessary to 
consider this matter further here.) 

Since the principal mass of the fluid is turbulent 
T^ and T^ are quite similar so that in general Tt can be 
replaced by the known 'Tu» -^t high Fr , tem.perature T^ 
agrees with T^ , hence with Tij {fig. 5); but at low Pr, 

is substantially loiver than T^. . Since is appli- 

cable only to the term of the transitional layer, the 
approximate value from; figure 5 will suffice. 

Of particular influence on the heat transfer is the 
temperature relation of the Prandtl number iii the laminar 
layer if a very viscous flu.id is involved. (Por viscous 
fluids the m.ajor resistance to heat transfer, is offered by 
the laminar layer. Since the laminar layer thickness ya 
varies v/ith the temp, viscosity history of the fluid (see 
equation 33) and since the temp.' variation of Pr is fixed 
primarily by the viscosity (cp and X .vary but slightly 
with temp.). Yg^ depends on Pr/Pi-Q. But it may not be 
said that the thickness of the lam.inar layer is a function 
of Pr for this'is a. heat transfer factor and the laminar 
layer thickness depends on a purely hydrodynamic variable 
as seen from equation 33^) In this event a m.ere estimate 
of Pr^ by means of a cursory tempe?a,ture T -j^ may intro- 
duce a seriotis error. The pred ic t i on' 6.f • * P.r by means of 
equation (o5) is t heref ore indicat ed. 

To simplify the calculation a proportionality betv/een 
velocity and temperature for the hon-^ i s o t he rmal lam.inar 
layer is assumed: 
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VdA 

* T = CP i^— (43) 

Pr 

Further the temperature variation of — is approximated 
by the linear equation: 



Pr Pro ■ • 

= (1 + m T) (44) 



v/here m is an empirical constant. 

With the assumptions (43) and (44), equation (35) gives 



0 



7 dt^> 

A first approximation to the temperature gradient ( 

at the wall is obtained from figure 4 where for. the 
Prandtl numVer Pr * the value of the wall, Pr'^ may 
he ch 0 s en • , 

?or a more accurate solution of Pr (to he discussed 
elsewhere) the real non-linear functions T ( '^) and" Pr (T) 
must be used instead of the linear rela-tions given by • 
■ equat ion^ (43) and (44) • ' ' 

■ THS DISTRIBUTIOIT OF THB HEAT ?LOW DENSITY 
IN A CHANNEL. AND A PIPE 

» 

'The differential equation for the equilibrium of 
forces . in a fluid with- allowance for the continuity equation 
and omission of density variations- reads:' 

p — + p V ('t7;w) ■ = \. grad p' -f. ' u w ! ■. (36) 

ri t ■ ■ • > '■ ■ • 
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where v is the velocity vector. 

The differential equation for thermal convection 
and conduction is v/ritten similarly: 

P T"^ + P Cp V (w;T) = e+ \ V T (47) 



v;here 

€ source of density per unit volume, 

A formal analogy hetv/een the equations of momentum 
and heat exist therefore for flov/s v/ith grad p ^ 0 
only in the presence of spatial heat sources in such 
flov;s. Even though the internal friction of a fluid 
is small in technical applications, the variable ^ is 
retained in the equ&.tion for future consideration of the 
analogy* 

Equat ions ( 45 ) and (47) are next applied to the com- 
pletely developed turhulent flow in a flat rectangular 
channel and in a pipeo Tor this type of flow the non- 
uniform terms cancel out hy averaging and likewise the 
derivatives of the mean velocity along the principal 
flow X, 

With u = u + u' and v* denoting the velocity and 
fluctuating compon.ents of the velocity in the x and y 
directions., respectively, where "bars represent mean vel- 
ocities v/ith time and the primes represent instantaneous 
variations from the mean, the scalar equation for two- 
dimensional channel flov/ in terms- of m.ean values is 

In this equation the "bars are omitted from the pressure 
and velocity terms as v;as done previously. The hars are 
used only to represent the mean products of fluctuating 
quant it ies . 



The moan momentum interchange — p u* v', which may he 
regarded as a stress attitude is identical with the tur— 
Dulcnt shear stress t^, while |a, is the viscous 



(48) 
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shear stress Ti. (t^ is the turhulent momentum 

transport in the — y direction. For ^ > 0 '^t is 

dy • 

likev;ise positive; Bince the positive' ar e ' as s o c iat ed 

v/ith the negative v' and vice versa on the average, 
^'or the exchange process the higher* u velocities arrive 
fro::i .greater and t ho ■ 1 6w v eloc i t ie s from closer v/all 
distances ,) 

Hence is the total shearing stress 

5 5p - , \ / \ 

— = — = constant (channel/ C^y; 

oy ox 



The pressure drop is constant since the flow is fully 
developed. 

2hc heat "balance' in tv/o— dimensional channel flow 
follows from cq^uation (47) 

\ bx by J by^ 



Assuming that no great changes in (T^ uO occur in the 
direction of flow x as certainly is the case for fully 
dcvolo'pcd temperature d istr ihut io-'ns the second term on the 

left side of equation (50) may he neglect odj^ Ther e rem.ains , 

then, only the fluctuation product, p Cp 5' v^ which is 
equivalent to the turbulent heat transport" perpen- 
dicular to the v/all» Introdvicing the total, heat flow 
ST 

q = Ox + X — , further affords: 
oy 

^ " P '^P ft ' (channel) (5l) 



Por. fully developed flov/ in -a pipe the . f o llov/ing 
equation in cyl indr i cal co or dinat e s -( ins t ead of equation 
(49) for the channel) is o"btained 



£1 =.|2 = constant (pipe)' (49a) 
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arxcl (inctead of eq^uation (5l) for the channel) the 
equationJf 

~ (q(l~T])). - (l^Ti) f P Cp n 1^ --i" (pipe) (51d) 
oy \ ^ ^^"^ . / 

where ri = ^ (r = radius of pipe or l/s v/idth of channel). 

3;/ introducing the ohear. stress at the wall 
there follows for the pipe and the channel frora equations 
(49a) and (19), respectively: 

T ~ ( I'-T, ) ( channel .and pipe) (52) 

Por the case of fully devel '^p od " t einpo r atur e distribution 
to v^hich the argument is confined the partial differential 
oqiiations of heat flow becone ordinary differential equations* 
(The entrance lengths for fluid flow and heat diffusion are 
independent of oa.ch other. Per instance, a sudden change 
of v/all tempera'cure maj introd-j.ce a thcriaal entrance in a 
hydr odynaLiically developed flow,) The s i:n ilar i t y of the 
t eiiiper a.tur e profiles states that the decrease of t emp or atuir e 
per unit t ei^perattir e — dl/T i?r. ■ di stance 6.x at each dis- 
tance y from the wall is a constant. Therefore — dT/dx - 
constcint T. The onission of e leaves • ins t ead of equations 
( 51) and ( 51a) 

-r^ = constant u T (channel) (53) 

ay 

d 

■ — (q(l--r, )) - constant ( 1— ri ) ti T (pipe) (53a) 
dy 



T is defined, as before, as the temperature excess over 
that' at the v/all. 

It should be noted, however, that a fully developed 
te:npcrature dis tr ib'ut ion is possible only at small excess 
temperatures T. In general u(S/c):c)(p c-p T) instead ef 

p c-Q u would have to be r eckon ed "w i th ^ Por the f^llow^ 

ing it is assui-ed tha.t change of profile remains v/ithin 



Z6 
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such limits that equations (53) and (53a) remain applicalDle 
v/ith sufficient accuracy. 

xZxZ e^T Ion o£ equations (53) and (55a) ana xu"cr oduct ion 
of the limit values ( q = qo at "n =0 and q = 0 at ri = 1 ) 
in v-^^on j unct ion v;ith the nondimensional ^ = T /9 and . 



= u/U r;ives; 



q. 



f <^ <tp 6. T] 

= 1 ^ (channel) (54) 

^0 / 1^ cp d ri 

. 0 



S ^ cp( 1-ri )d r\ 
^ (l~Ti) = 1 - '^^ (pipe)(54a) 



/ ^ (p( l~r, )d r\ 



o 



A t-^ood appr oxiiTia t ion for q/qo ootainod v;ith the 

a2:'plico.t ion of the temperature equations (37a) to (37c). 
Although the introduction of the simplified temperature 
equations (40a) to (4l) is sufficient. 

Tho velocity 'distr i out ion of the turbulent region 
is represented "by the v;ell knoi^in pov;er lav; J 

9 (55) 
v/here 0.18 > ri > 0.10 for 4.10^ < Re < 4.10^. 



The pov/er lav/ represents the velpcity distribution 
even better for large T) ' s than the logarithmic law and 
is especially suitable for the present calculation. . In 
•wall proximity the velocity is of course less than that 
calculated by the power lav; and the error of the derivation 
is small only at high Reynolds numbers where the power law 
must be used near the v/all. 

The ur. 0 of equations (40a) to (4l) -and equation (55) 
then yields appr oi-cimat ely : 

, (1 + n- n2)(Pr» - 1)9,^+ 9 a+i/n . 

/ 4 9 d ^ -^-7 — ; ^ —sJP^ (56) 

J (1+ 2n)((Pr' - l)cpa"*" 1+ 0) 
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.Here the lower limit cp^^ is replaced hy 0* Because the 
value of this integral is practically r^ero at the, limit cp^ 
since the ..exponent .1 .+ l/n is high and the integration 
from 0. to cp^ y i elds .an - int egr al ///h ich is very small* 

2hus the density of heat flow for channel :f low -at 
high Reynolds numbers approximates toS 

g • •• (Pr»-.l)(l+n)cp +cp .^^/r> , w n 

— = 1 - ■ ^S:—— cpi-^i/ ^- ( channel)( 57) 

^0 . (Pr »-l)( H-n)cpg^+l 

n , ■ 

v/nore cp- = r) , so that q./qo represented as a 

function of . ri . . 

In a similar fashion the heat flow through the fluid 
in a pipe is approximately! a.ccording to equation (54a): 

— ( l-Tl ) = 1 

(Pr»-1)( l+2n)( 2+n-( l+n )cp ^ )cpa+( l+n/2)( 2+2n-( l+2n )cp^'^)cp / 
: . __ .— cp ^+ ^ ' >^ 

(Fr'-l)( l+2n)cpa+l+n/2 

(pipe) (58) 

A good view of the variation of q/q.0 ^ay te obtained 
for the special cases of Pr ' = 1 and Pr* — > co as sud~ 
stitutod into equations (57) and (58)': 

■ = 1 - Ti • (channel Pr' = l) (57a) 

q i + n , '. ' V . , V 

~ = 1 - Ti ■ (channel Pr' = °°) ( 57h) 



1. (i^n) = 1 - (2 + 2n - (1+ 2 n)n)rt'''^^ (pipe Pr« = l) (58a) 

"0 



q 



(1+ri) =.1.~ (2 + n - (1 + n)T\)r] (pipe Pr» =00) (58h) 



0 
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The o/qq curves flatten out with increaf?ing values 
of Pr and ?we (decrea?iing n), that i pi ^ they approach 
q/o^ = 1-11 which is the limit of equations (57a) and 
(53a) for 'n = 0. It -airther follows that the quantities 
k and e decrease "'•^ith increasing Pr' and He (see 
equations (13) and (39) and fig. 10). 

Figures 7 and 8 show (q/a^) ( cp) and (qJoq) (t)), 

respectively, for the pipe and channel at Re = 4 x 10"^, 
Figure 9 reveals (q/q )('n) for various .values of Re 

0 

at Pr^ = 0^72 and Pr ' = 200. These curves were com-- 
putcd for "by the true velocity distribution as shown 

in figure 1 rather than power law. 

The variation of the heat flow density In prOTimity 
of the wall is not eworthy, where (da/dT])^ = 0 for the 
channel, hut (dq/dri)^ = q^ for the pipe (see equations 

(53) and (53a)). The rise of the heat flow density of the 
pipe heyond the value q^ is due to the fact that the 
total heat flow Q, ^ qd-ri) near the wall is practically 
constant, while the section through which the heat flows, 
decreases with (l-n). In channel flow no cross-sectional 
area changes occurj thus the h^-at fluY density and the 
total heat flow are always directly proportional . 

At mid-channel .(and pipo, respectively) the variation 
Of the heat flow density is characterized by, 

/da^ ^ gp 
vd-ny^ -oucpm 

' '^u "^m 
(see equations (54) and (54a)), -^.^ = — and = 

denoting the d i men s i onl e s magnitudes of the mean stream 
temperatures T^ and the mean velocity ("^ir,)! respectively. 
Since these magnitudes are smaller than unity, the negative 
?^lope of q/q^ is greater than unity, and is greater for 
the pipe than for the channel. 

Equations (57) and (58a) enable the calculation of 
k and ei (see equations (13) and (39)) through which the 
second appr.ox imat i on to the temperature is secured. In fig- 
ure 10 the ei term for pipe flovr is shown plotted against 
Pr« for different Re. 

(As mentioned above (equation (9)) I. Taylor has 

computed a second approximation to the temp, distribution 
for the case of Pr = 1. His arguments rest on Reynolds' 
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analp^-^y. .In the calculation dT/3x ^ is as.sumed inde- 
pendent of the dii=tance from the wall which' corresponds 
to, Bq/9y proportional to u (instead of u T). 'Through- 
out, a correction term, which is too small,' results (for 
instance at He = 4^10'^ it is 5 percent instead of 9 
percent, as is the case for high Pr" (that .is, for T ^ 
constant)).. (See fif. 10.) ' ' " ' ^ ' 

COMPARISON BETWEEN MOMENTUM INTERCHANGE 
AND* HEAT- TRANSFER 



Supplemental to this theory an attempt is m*ade to 
compare the differential, e.quations of heat and momientum 
and to indicate that t'he historic heat source theorem 
also leads to a generalization of thePrandtl humter. 

The sim.ilarity of the differential equations (46) 
and (47)" is so "obvious tliat 'if need not be discussed - 
further. However, it is necessary to analyze the 
existing differences, 

One sub st ant ial departure lies- in l^he fact that t'-he 
heat -equation contains no term corres-ponding to the ' 
pressure drop in the m.om.entum equation.-- This difference 
can, however, be removed in some Cases (as Prandtl has 
shov;n) in first approximation by substituting a suitably 
'chosen heat source density ^. 

Physically this artifice has the following signif- 
icance. ..The m'om.entum of a flowing fluid can be maintained 
by a pressure gradient. The hf--at content of the fluid 
i s, in • c o'nt r as t , reduced by t-he transfer through. the walls, 
unless heat is produced in thr fluid itself (such as by 
a current of el f^ctr ic it.y '.flowing through thefluid). To 
c-omr;lete the analogy between .heat transfer and m.omentum 
exchange the volume -heat sources must be s o di spos ed " that 
the t emperatur e and velocity -.prof iles are ; s imi 1 ar , In the 
particular Case wh^-re the velocity d i ? t r i but i c^n .•■ r ena i n s 
constant in the axial direction (fully devel oprl flow in 
a pipe or channel) the temperature profil-F' should be main-* 
tained likewis'e, 

A further . di ff erence between equations (46) and (4?) 
rests on th^-fact t hat e quat i on- (46) is a vector equation 
and equation (4.7) i.s a scalar equation, hence only one 



I 
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component, of eouation '(46) can be, compared with equation 
(4-7). Th'e a'n'alogy i s , therefore,., carried out for fully 
develop^ed plane .channel; ^low-. "By. concentrating on the 
sp'eclal case. oT. plane; frow,- the problem is much clearer, 
and affords mor-e far reaching co*nclusions than from 
equations '(46) and .(47)., 

For this flov.^ equations (46) and (47) "give: 
d / ■ du — ; -N dp 



(49b) 



^ S dT — -\ ^ _ V 

i7 " V,^T, - f"=p ^' Vy = - c . P s u - C61.) 

In order to avoid misinterpretations', the mean values 
are again represented by bars. 

In the most freneral form, tiie equations for- T 
and a are analogous, But because fully developed 
flow has 'been postulated, the terms vfith du^/dX', 

5 -i"- • ' • • ■ . : . S 

dy^ (T. -v)" ■ drop ou*t , while the term with (uT). 

in the heat equation remains. The p.ro'blem, then, is to 
choose c so that the term "with ^Y/'dx vanishes* 

As shown above • dq/dy*^ -u St/Sx,^ u T (see equation 
■,(5?)). _ Thus for a given heat. source d i s t ri'^iit i on_ 
£ = k u T, a suitable choice of k 'will cause dT/dx 
to diminish to zero- for every y, without^ in any way 
modi f y ing the temperature distribution T or the exchange 
T^ v' , With this choice of' c the' actual temperature 
profile v;ill be retained; it 'simply r^emedies the earlier 
decrease in tem.perature in the direction of flow. 

However this heat density -which varies with the 
distance from the v/all cannot be compared with the 
pressure gradient which is constant over the section. 
So in. order to carry through the consideration of the 
analogy it is necessary that ^ - const_antr In this 
caf^e also it may be stipulated that dT/dx should 
vanish at each position of y. Then the heat flow is2 
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— = - G = constant 

Under this condition T is of course no longer the 
actual temperature, but rather an approximation to the 
temperature, which is that the postulate c = constant 
yields too great a. temperature gradient at the wall he- 
cause the'V/all layer? are heated exc-essively "by constant 
source density. But in view of the fairly well compensated 
temperature profiles the error must he sma.ll. 

iTow the identity of epilations (49"b) and (Slh) can 
"be adduced hy putting conformably to Prandtl, [j, = X 

•(Pr.= l) and T = 3 u^(p = constant). 

Put the equation T - p u is only one possible 
solution. The solution is, in fact, somev/hat special, 
since it not only requires the time averages of the 
velocity and the temperature proportional, but the 
fluctuations u' and v' themselves to be proportional 
at every instant. 

T = P u • • ■ ■ . (59) 
T» = P u» • (59a) 



These equations are obviously fulfilled if the mech- 
anism of transfer of the u-component of the m>omentum and 
the mechanism of heat transfer rt^q c-ompletely similar, 
. This .may occur in particular oases.' 

"Consider n^j^t the general case where the correlation 
coefficient between u and T is less than 1. To this 
purpose the turbulent terms in equations (49b) and (51b) 
are expressed by .A du/dy and cp Aq(dT/dy), The 

identity of the equations is attained vrhen A - Ai--2^ 

(that is, Pr » = l)' for- each - di s tance y from the wall 
and when equation ( 59 ) . i s s at i s f ie d . 

Equation (59) thus re;presents an approximate solution 
of the heat equation for the case defined 'by Pr^ = 1» 
In .this case a.^sumptions relative to the fluctuations of 
u' and T' are no longer required. 
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The solution (cqu^aion (59)) stipulates that the fic- 
titious source strenirth € b e f ix ed "by : 

■ <^°> 

Cn the basis of the postulate dq/dy = - ^, the 
source strength € may he defined in terms of heat 
flow at the wall q.^ : 

e = q^/r (61) 

(r is the half*-width of the channel)* Introducing the 
shear stress at the wall for the pressure drop, it 
follows from equation (60) that: 

Aq 

^0 = P ^-P T ""o . (^2) 



(This change of form 
that the form of the 
plate, and so forth) 
siderations does not 
may he derived also, 
pipe. The constant 
standard values (for 
the miyed mean tem:-p. 



of thp equations has the advantage 
fluid "boundary (whether pipe, channel ^ 
which is unessential for these con- 
affect the result.) This equation 
for the qua si "-plane case of the 
|3 can "be expressed "by appropriate 
instance, hy the maximum values) 
of T or a,. 



At Irr» 
total fluid, 
the terms of 



/ 1, the analogy is not comrclete for the 
"but only for the turbulent- core in which 
molecular conduction can he neglected. 



Then equation (59) is appropriately replaced hy 



T ^ Ta = & (u - Ua) (59h) 

where T^ and u^ are time averages at .^Hhe point of 
transition to the laminar flow".- For the rest the 
Calculation is the sam,e and equation (62) holds for 
■Pr ^ 1. "Solely 3 becomes another proportionality 
factor. 

Jor the prediction of p the maximum values of 9 
Or U, the temperature and velocity, respectively, can 
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iDe \ised. It is also posoilDle, however, to introduce the 
sectional averages Tj^^ and Uj^, (as was done in the 
Prandtl derivation) since the section of the laminar layer 
of the flov; is negli£:ible compared v/ith the total flov/ 
section. With T^^ expressed in terms of Ug^ and Pr ^ : 

T^ = P Pr» Ua (63) 



- an equation 7/hich results upon the application of equation 
(62) to the laminar layer, equation (59"b) gives: 



or 



3 = ^ (64) 

1 + (Pr »-l)u^/U 



p = ' ^ (64a) 

1+ (Pr »-l)ua/u^ 



The following useful conclusion can "be drawn from 
equations (64) and (54a): 

co^ + CD (Pr»-l) 

^m = ^~ T (65) 

1 + cpo (Pr »-l) 



T u 
m , cp _ m 

0 ' ~ U 



The dimens ionles s mean temperature is 



identical at Pr ' - 1 ' v:ith. the d imens ionle s s mean velocity 
(v;ithin the framework of the present approx,), V/ith increased 
Prandtl numher approaches unity. 

Since it has been established that the correct heat 
flow distribution over the section is not necessary for 
the determination or .the temperature distribution in first 
approximation, q/T = q^/ and this assumption compared with 
the assumption of a constant source strength, which as 
shov/n above, is necessary to establish the analogy. From 
dq/dn = - r = constant, it follows that q = q (l-ri)* On 
the other hand, for completely developed flow t - (i-p,)^ 
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For fully developed flow our approximations therefore 
agree with the postulate of constant heat source distrib- 
ution. The two idealizations differ, however, in their 
consequences for further theoretical treatment as well 
as in the justification of their physical admissibility, 
although "both methods of treatment agree with great com- 
pleteness as a "basis of the turbulent velocity profile. 

(Since the heat source theorem is designed to describe 
analogous phenomena which in detail are not analogous, the 
theory is. definitely bounded, which limits its extension. 
The heat source theory is not purposed to consider particular 
force fields in detail, merely intended to reproduce and 
clarify the essential characteristics of the heat flov/ by 
compp.rison with known phenomena of the momentum transfer.) 



HEAT VOLUME THANSMITTED TO THE WALL 

(a) Determination of the heat transfer from the 
temperature gradients at the v;all. 

From the temperature rise at the wall the unit heat 
rate at the wall follov/s directly at: 



This equat ion , -1. : naturally general, is not affected by 
the type of flow nor by the time or space variation of 
Tq and (du/dT)^, 

The rate of heat transfer Q, over the area ? is 
obtained by integrating equation (55). If U is the 
velocity and 9 the temperature difference of the fluid 
v/ith respect to the v/all at the boundary of the friction 
layer, the heat flow is: 




(67) 
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Por p.n area over v/hich • 9, U, and (d^^/d^^)^ are 
sensil)ly constant it follovjs further that: 



^ = Jt2o_:_ ^ (68) 
^Vd-Vo 

where 

AV/ rer.istance to flow offered "by the area under consider- 
ation. 

If 0. turbulent friction layer is involved Pro (d^^^d-^)^ 
in equation (68) is expressed by equation (SS). 

As to the permissihle size of the area to which equa- 
tion (68) can be applied in friction layers free to extend 
unhindered over the surface (that is, the actual "boundary 
layers), © is, in ^.eneral, the constant temperature 
corresponding to potential flow. In this c^use the ad- 
missible size of the area is dependent on the adequate con- 
stancy of U and (dcp/ d.^) ^ CD^ and co-^, respectively, 

of equation (38) - cp^^ a.nd cp-^ are functions of the Re 

of the velocity profile and so vary v/ith the arc-length x) . 

In flow through pipes or channels, the Reynolds number' 
for fully developed veloc ity profile is constant; but the 
maximum temperature © decreases, 5^or which reason equa- 
tion (68) holds only for pipes and channels if the flow 
section is sufficiently -short* In long pipes the tempera- 
ture drop must be accounted for, as shown in the next sec- 
tion* 

(b) Consideration of the heat loss in the friction 
layer. The heat stream Q,^ flowing through a flow section 
f is 



Cp T u df 

Introducing the mean flow temperature (that is, the m.ean 

temp* of the fluid mass flowing through the flov/ s ec t ion) affords 
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f 

/ T u d f 



'0? = 



u 1 
m 



or for constant values of Cp and p: 

" P ^ 

Che heat volume given off at the v/all over the arc 
length (x— x^^), is equal to the difference Q = Q,^ ~ ^ 

in the ilov/ sections f and f^. V^ith the postulate 
that the profiles of the velocity and t e-^-por atur e are 
siiuilar: = .^u/®» ^^'^ ."^mj'^i = ^m/^^ hence 

•V Q = P ^P --a 9ri (Ui ©1 fi U e f ) (59) 

v;horo ^ T^/S and 9n = ^^n/^* 

Equation (59) can "bo utilized to check experimentally 
, the theory (similarly comhined v.'ith equations (38) and 
(5S)), All quantities in equation (69) are readily measure— 
ahlo, A minor complication is introduced in fixing 
for vfhibh a -mixing chamber is required. • 

Also ^ can be evaluated from, the theory^ Per the 
plane case: 

/ : 1 

o '"^ Cp d T] ^/ cp d Ti (three dimensional 
= friction layers) 



^ 9 d Tl 



and for the flow in a pipe 



1 1 

/ . ( ) d Ti / cp ( 1-Ti ) d ri 

0 V 0 

=: 2 ' = — ; ~ (pipe) 



0 



In figure 11, '^^ is shown for the flat plate and 
for the pipe as a function of Pr for several Reynolds . 
numbers^ - Tho effect of the PwCynolds number on is 
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less" than the effect of Pr , since cp(Rg) appears in 
the numerator and the denoninator* 

(l) Pipe or channel ilowr 

In flow through a pipe or chcannel, the development 
of the friction layer is limited and the strength of the 
final friction layer is equal to one hraJ.f th distance to 
thQ opposite wall. The mass of fl'uid moved m the fully 
developed friction layer does not change, and the heat 
given off ty the fluid can "be calculated from the reduction 
of temperature^ 

So, when no change in cr o s s- s e c t ional area is con- 
sidered equation (69) gives; 



p Cp ^.^ u^ .Oi-B)f (70) 



On the other h n d , t h e . h e a t transferred is also de- 
fined "by equation (67). Here, it must he noted that, he- 
cause of the similarity of the temperature profiles the 
percentage temperature drop dT/T over the len.?:th dx 
is the sr.me at all distances from the v/all^ Hence it 
also applies to the mean temperature: 

In — = constant (x-xi) 
and equation (67) yields 



In [—J Prol—j 



v/he r e 

given ma^ximum initial tomper.ature 

V/ frictional resistance of the pertinent pipe or 
channel length. 
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(At sniall tempor.^turc differences 9i-6 equation 
(71) changes to equation (68)). 

Sy equating equat i ons ( 70 ) and (7l) , the reduction 
of the mean temperature is: 



rlere the d imensionles s variable a* defined in equation 
(82) is introduced and s is the perimeter of the section 
f . 

?or pipes s/f = 4/d for flat rectangular channels 
s/f = 2/h (d = pipe diameter, "b = channel width, h = 
channel height, b»h), 

(2) Boundary layer flov/. 

If- th.e flow along a wall is not "bounded, the friction 
layer can develop unhindered, and v/hile the boundary layer 
increases in thickness in the direction of flow, the maxi- 
mum temperature 9 on the boundary of the friction layer 
remains^ in general, unchanged. 

If the surface 01 the body has the temperature 6 of 
the fluidi then the heat flow density at a particular point 
in the friction layer is equal to p Cp 0-^ . But the cool- 
ing auction of the v/all lowers the temperature by B- T, Th 
"cold -.stream" through the section f of the boundary layer 
at the point in question. ..is therefore: 



In 




(72) 



f 



ft = 




u d f 



(74) 



0 



where 



heat volume absorbed by the body surface up to the 
particular point x in unit time. 



^y introducing the dimensi onles s value of the mean 
flow temperature, equation (74) gives.: 
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q ="p" 9' (l-x^^)f (75) 



the value of -e^-^ for any jR;iven velocity profile is 
approximately defined "by the theory. 

At constant maximum velocity U the heat alDsorlied 
over a length x-x, is. fixed "by the increase in section 
of the boundary layer thickness in the direction of flov/,' 
since cOjj^ and vtiry but little {^^^ increases and 

(l-i^-u^) decreases with the Re of the boundary 'lay er ) . 

The effect 'of Pr on the heat transfer is expressed 
by the factor (l-^^). With increasing Pr " the tempera- 
ture profile becomes more blunt-nosed and -o- approaches 
unity. 

Prom the momentum lost in the boundary layer relative 
to potential flow the flow resistance W of the body can 
be written in a similar manner to that used for the heat 
diffusion Q, 



• \I = p U u,^ (l-^\)f - (76) 

f ■ 

Eeroin cp^ is defined by cp^ ^m - = d f . The differ- 

ence between the equations lies in the velocity U in • 
potential flov/ which, in general, is not constant like © 
but v?.ries v/ith the arc-length x. . . 

But the analogy betvj^een equations (75) and (-76) is to 

be carried out f-or a surface area over which .6 and U 

are constant, (a pressure' drop is to be avoided). If . 
f-f^ is elimi-nated, and if and are the mean 

values over the particular arc length; 



= Cp 3 i — 77 

Here the heat transfer is expressed by the heat loss in the 
fraction layer, while in equation (68) it appears in terms 
of the temperature gradient at the wall. Comparing equations 
(68) and (77) yields: 
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1 - CP 

\ = 1 (78) 

Jot- = Pr^ ^ I, (dcp/d^)^ = 1 and hence = ^D^, 

which is due to the similarity of the velocity and temper- 
ature profiles for ^ero pressure drop. P'r^ (dcp/di-OQ in- 

cr.e'ases with increased Pr and ^« 

■ Por the rest, equat ion :(78) is. easily verified for the 
simple friction layer (equations (40a) to (4l)). 

(c) Keat transfe.r coefficients. 

The heat transfer coefficient a is defined "by 
Newton's law of cooling: 



= a (T-T.) AP (79) 



in v/hich the heat transfer per unit time through the 
"boundary area A? is put proportional to the temperature 
differe-nce. (T.- .Tq). of the fluid and of the wall* 

Originally the proportionality factor a was thought 
of as a pure m.aterial value comparahle with the thermal 
conductivity and in the older literature v/as designated as 
the "outer thermal conductivity". V/ith the increase of 
experimental data, it "became more and more apparent that 
the flow phenomena ad jacent • to the wall contri'buted greatly 
to the 'heat transfer and varied in a complicated manner 
therewith. Hence the cooling law is only apparently simple, 
that is, when the. simple form of equation (29) is maintained, 
all of the pro"blems of heat transfer "by convection are con- 
densed in the factor, a. 

This naturally does not. help to clarify the physical 
phenomena and later re-search has produced other axioms which 
throw light on the mechanism of heat transfer. For the prac- 
tical" application of research data it is, however, advantageous 
if complicated relations can ho expressed hy a single coeffi- 
cient^ the value of which is ohtaina'ble from graphs or tables* 
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In equation (79) the fluid temperature T remains 
undefined.. What temperature between Tq and ^max 

use is purely a matter of expedie.ncy. Only one point is 
necessary, namely, that the fluid temperature employed is 
adequately defined. The reference temperature m.ust "be 
.relatively constant over the area under consideration.. 
Since this condition holds true in all cases only over a 
small area AF. the heat transfer coefficient must he 
defined as local quantity.. 

In the earlier derivations the heat transfer coeffi- 
cient was "based on the mean temperature 

The calculation of a heat transfer coefficient 
^m ~ ^o/^m means of equation (66) requires an express- 

ion' for — -L^^ from the theory, if is referred to 

(dt)^ 

the mean velocity u^^^.. Such an equation can he obtained 
from, equation (37c) after forming- a mean value of the 
velocity and' the temperature over the section of the 
turbulent region. Since this section is not much smaller 
than the total flow section it can be approximated to: 

Pr^ ^ ^1-H e^^^ i^,.^i)^^^^(llLlll . i\ (38, 

VdT/o ^m U - ' 

\ Fr 



This equation applies so much more as the boundary layer 
is thinner, that is, as the Reynolds number is greater. 

Neglecting the term. e^, and putting u^ = u-^^, and 

- A, equat i ons ( 66) and (38a) give an expression for 
^m ~ ^o/^m which is identical with the Prandtl formula 
e quat i on Tl) 

The mixed mean "temperature T^ is usually employed 
rather than the mean temperature Tjn v/hich is difficult 
to measure. This requires an equation for = q^/T^ 

from the theory. Here a difficulty arises. Proper treat- 
ment of equation ( 37 c )• af f o rd s a formula for (du/dT)^ 

in which the mean temperature of the turbulent region is 
the reference temperature. Eut then the equation includes 
a mean square value of u over the turbulent section in 
place of the mean velocity across the section. 
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By ap-pr oxirnat ion the r.ean of the squares of. vl 
can, 01 course, Id e replaced "by m^^^ (particularly at 

high Re) just as - -T^ can he roughly approximated in-- 
stead of But this also neans returning- to equa- 

tion (58a), that is, the Prandtl formula. 

The difference hetv/een and Tjvj played no part, 

hov/evor, in the earlier cons ider at ions, . ^ In view o'f the 
experimental difficulties, this difference usually lies 
v;g11 v/ithin the., experiment al . error . Further, the omis- 
sion of the transitional layer and the postulate of in- 
variable material values accounted for larger discrep- 
ancies- than this temperature difference. 

In developing a theoretical equation for qo/^u 
v;hich is in accord v/ith experience, ten Bosch (refer- 
ence (l?)) proposed a semi— empir ical equation^ The form 
of the Prandtl equation v/as follov/ed, hut the constants 
v/ere replaced "by variables v/ho s e* magnitud es v/ere deter- 
mined as a function of Re and Pr from -measurements 
available, the resulting equation "being* 

qo 0.125 i 



(for heating P 1.4 for cooling p ^ 1.12, Pr^ refers 
to the layers noar the wall. In the remainder of the form- 
ula the -properties are fixed at th c mean f 1 ovr t emp . ) 

In connection with the theory of the present roport, 
the heat transfer coefficient is most appropriately expressed 
in terms of the maxiinuin temperature 6, since only this 
temporatiire can ho used as a reference quantity without 
r ose:.^vat ion (see the derivation of equation (o8))» The 
maximum temperature ho.s in addition the advantage that 
it can "be measured withaut the use of a mixing cup. (In 
cases whore can ho measured m.ore reliably than ©, 

6 maybe calculated as ^u/'^u*^ Hence th o ^dof in it i on 



(80) 



A Q 
9 L r 



e 
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or in nondimens ional representation "by dividing "by 
Po^po^m* • 



a, = ^ ^_ ^33^) 



Po <^po ^ ^ ^ Po °po ® 



In the technical literature the heat flow is 
usually referred to p Cp u^^ (see equat ion ( 80) ) . 

Further, it is customary to Introduce the Nusselt numher 

llu = d/X 

in v/hich q^ is referred t'o as the heat flOv\f A -^u/^' 

Those two dimens ionless heat factors are related as 
foil 0¥s : 

^o Nu i\Fu a* 



where 



p Cp Uj^ He Pr Pe ^u 



p c-p u d 
Pe = L ^ ni = Peclet numher 



^ The older Nusselt numher had proved itself in the 

representation of cases v/here the heat transfer phenomena 

were not to he separated. But, for those cases where 

statements can he made relative to the local heat transfer, 

and where q^ can he written directly in terms of 

A (dT/dy)^ it is proper to refer the known quantity q^ 

to the product p Cp u T as will "be seen in suhsequent 
derivations. 

In certain cases it may he desirahle to compare the 
two dimensionless groups* 3ut»in generali both Nu and Pe 
are superfluous if a well founded formulafor the dimension- 
less group qo/p ^p ^ availahle. 
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In order that-" Q^/o\ 9 need not "be repeated 

^0 ' 0 po m 

unnecessarily the symbol a* has "been intro'duced. The 

dimensionless factor a* refers the heat volume 

transferred to the wall to a heat volume Cp u^ 3 

flov/lng past the wall. The njmlDer a* is, therefore, a 
kind of "efficiency" of the heat transfer. Ordinarily 
a* is very small (10-^ > a* > lO""^, see fig. 13) so 
that only a, small fraction of the hep.t becomes useful 
for transfer. 

^ Since a* represents a locally defined quantity 
as well, equation (68) Tnust he- "applied for the subsequent 
treatment of equation (8l). Then equation (19) in con- 
junction with cpjjj = "^m/^^ give 



(82) 

8 Pr (^) 



This formula,, is a.s general as equation (68) and the 
definition equation (19),- It is the'refore applicable, 
independent of the character of the flow, 

Pormula (82) had already been utilized to introduce 
a* in equation (72) in order to establish the tempera- 
ture drop which accompanies pressure drop. On the other 
hand, equation (72). can eq_iially be used to define a*. 



v/hich for a sufficiently short section affords 



This equation is, with regard to equation (70) identical 
with the definition (equation (8l)) for a* , 

Jot turbulent' friction layers the Pr^ ( d cp/ d -^^) q*' from 
equation (38) must be introduced in equation (82). With 
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''■-^''7» ^^^^ = - ~ -and assumin-g Cp^ - ^-pu ^pt 

U • ■ 

(which is permi s s itle ) 



0.1-5 i ep^ A^/A 



a- = (85) 



UVA • W • U \l-A/Aq_ Pr- / 



v.rhere 

t> m » U/u* are knov/n functions of He 

In equation (83) all variaoles and coofiicionts are 
known except a.-D," and -^n/-^^' flow- measurement s 

a = 1.5, "b = 15.5. The. vdlue of 1^55 may' l^e used since 
"b. occurs only in' the term referring to the transitional 
layers But the assigned value of a is very uncertain. 
Hence a and A^/A must oe determined from measured 
values of a?* accordin^-^ to equation (33). • 

Admittedly there exists a certain difficulty, involving 
two unknowns but equation (80) indicates that Aq/A scaresly 
■affects a* at high Pr, and a has "but little influence 
on a* at low Pr. Hence a can "be from heat transfer 
measurements at high Prandtl numbers and " Aq/ A from 
similar m,easur ement s at small Prandtl numhers.- 

Pi^-.-ure 1? illustraterr' the measurements by 3uhne 
(reference (18)) and Morris end Whitman (reference (19)) 
ait high Prandtl num.bors and the measurements by Hohonczi 
(i.e.) at low Prandtl numbers compared v/ith predicted ex*. 
Por the calculation a = 3.2 and A^/A = 1 were chosen 
to achieve the best correlation of predicted with exper- 
imental results. The spread of the test points at high 
Pr is understandable because of the difficulty of measure- 
ment at high Pr and the fact that non-isothermal flow 
theory is not yot complete. 

Por constant or slightly variable material values only 
one value of Pr enters in equation (83). Por this case 
(fig. 13) shows a* as a function of Pr * at several > . 
Heynolds numbers. With the help of a* established by 
equation (83), the heat diix-asion LI of a given area 
A P c an be c a. 1 c u 1 a ted. 
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= Po ^po % 0 A 5^ a* (68a) 

This is a form derived from equation (8l) v/hich is 
another form of equation (68).- 

The variable a* is a point f-unction and holds only- 
over a small area,' The area must he chosen' Just large 
enough go that 6 is sensibly constant over it. 

Equation (68a) .cannot, "be used for long pipes in 
which the temperature changes materially. Here the heat 
transfer may be computed by means of equa-tion (70) if 
the temperature reduction (8i-6) in the pipe length 
under consideration is known. The end temperature Q 
can be predicted from equation (72) when the initial 
temperature ©i is given. 

As is seen, the variable -^^^ likewise plays a relet 

But- this is only true for long pipe lengths for here the 
heat loss in, the_ friction layer is decisive. Por short 
pipe sections does not enter (this is also the region 

in which the conductance is applicable) (see equat i on ( 68a) ) • 
Then the t em.perature s are practically constant a.long the 
short, surface length, and the heat diffusion A(i is given 
by equation (68a) as a function of the temperature gradients 
at the v/all, v/hich fact is included in egua.tion (83) for 
a* 



Translation by L. H. K. Boelter, 
University of California, 
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greatest closeness 
are not sufficient 



"surface tube;" probably attained 
to the wall. But even- these test 
for the present con s ide rat i on s • 



the 
data 
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Fig. 1 
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Figs. 2,3,11 
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Figs. 4,5,6 
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Figure 4,- (d^/d^)p, 
the ratio 
of the temperature 
gradient (d^/dy)o to 
the velocity gradient 
(d<j)/dy)o both measured 
at the wall as a 
function of the 
generalized Prsindtl 
numbei- 




0 A3li Qi 



Figure 5.- Temperature (^) as a 
function of velocity 
{op) for flow through pipes and 
along plates with Pr' as the 
paraimeter. Re = 4 x 10^. 



Figure 6.- Temperature {^) as a 
function of ti = (y/r) 
where y = the distance from the 
wall for the pipe with Pr' as the 
parameter and the velocity dis- 
tribution ^ s ^ at a flat plate 
for Pr = 1. Re = 4 x 10"* 
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Figs. 7,8,9,10 




to 

Figure 7.- Heat rate q/qo and shear 

stress (r/ro) distribution in a 
pipe auad in a channel as functions 
of 'the velocity cp , witn Pr' as the 
parameter. Re = 4 x 10^. 





Figure 8.- Distribution of the heat 

rate (q/qo) and shear stress 
(T/ro) as a function of the wall 
distance rj for Re = 4 x 10"^. 
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Figure 10.- The variable p, as a 

function of Pr with Re 
as the parameter. Flow in a pip© 
with a = 2.0. 



Figure 9.- Distribution of the heat 
rate q/qo as a function of t) in a 
pipe for Pr' = 0.72 and Pr' = 200 
with Re as the parameter. 
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Figs. 12,13 
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Figure 12.- Heat transfer coefficient a* in tubes from the 

measurements of Biihne, Morris and Whitmsin and 
Rohonczi compared with the predicted values from equation 
83. (a = 2.2, Aq/A = l.l) . 
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Figure 13.- Heat trsuisfer coefficient a* as a function of Pr 

from the theory with Reynolds' number as the 
parameter (properties invariable, a = 2.0, An/ A = 1). Prandtl»s 
predicted results are included for a « 2,0 (dashed) and a = 8.8 
(dotted) . 



